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ABSTRACT 

The  dynamics  of  a  ship  propulsion  plant  are  modeled  in  an  all 
digital  simulation.   This  model  is  combined  with  that  of  a  mariner  hull. 
The  behavior  of  the  speed  governor  is  studied,  and  an  external  feedback 
loop  is  added  to  provide  direct  control  of  ship  speed. 
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I.   INTRODUCTION 

The  propulsion  system  of  a  ship  is  a  complex,  nonlinear  system  with 
several  internal  feedback  loops.   The  performance  of  the  ship  itself  is 
dependent,  of  course,  on  the  capabilities  and  limitations  of  the  pro- 
pulsion plant,  which  may  limit  overall  ship  performance  so  that  the 
capabilities  of  the  hull  design  are  not  realized. 

The  first  objective  of  this  thesis  is  to  model,  in  detail,  a 
typical  ship  propulsion  system,  and  to  simulate  it  in  the  digital 
computer.    After  validation  of  the  computer  model  a  simulation  of  the 
Mariner  hull  is  to  be  coupled  to  the  propulsion  plant  so  that  control 
system  studies  may  be  undertaken.   As  part  of  the  validation  procedure, 
behavior  of  the  shaft  speed  was  investigated  with  and  without  the 
governor. 

The  second  objective  was  to  use  the  computer  model  as  part  of  a  ship 
feedback  control  system.   For  this  purpose  it  was  decided  to  explore  the 
use  of  a  feedback  loop  to  maintain  constant  ship  speed  in  turns. 


IBM  360/67  DSL  360  LANGUAGE  at  W.  R.  Church  Computer  Center,  Naval 
Postgraduate  School. 


12 


II.   PROPULSION  PLANT  DYNAMICS 

A.   SHIP  PROPULSION  EQUATIONS 

The  propeller  angular  speed  is  denoted  by  n. 

The  basic  equation  describing  the  angular  acceleration  of  a  propul- 
sion shafting  system  is: 
EQ  =  2-irIdn/dt 

Where:   EQ  =  Summation  of  torque  =  Qe-Qp 
I  =  Polar  mass  moment  of  inertia 
n  =  Propeller  angular  speed 
Qe  =  Shafting  torque 
Qp  =  Propeller  reaction  torque 


then    dn/dt  =  -r~z   EQ 
2ttI   x 


For  the  typical  gas  turbine,  the  experimental  data  on  the  power 
turbine  torque  Qe  versus  the  power  turbine  speed  is  shown  in  Fig.  Al. 
The  torque  map  of  Fig.  Al  allows  the  engine  torque  to  be  determined  if 
the  engine  speed  and  fuel  flow  rate  Wf  are  know.   However,  this  torque 
representation  is  correct  only  for  steady  state  contitions. 

Propeller  speed  n(rps)  is  related  to  the  engine  speed  N3(rpm) 
through  the  reduction  gear  ratio  Kg 

N3  =  60Kg.n 

The  block  diagram  for  the  ship  propulsion  system  is  given  in  Fig.  A2 
The  power  turbine  speed  can  be  controlled  by  a  governor  acting  on 
the  gas  flow.   The  changes  in  fuel  flow  and  speed  that  occur  can  be 
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assumed  small  enough  to  permit  the  use  of  constant  coefficients  in  the 
dynamic  equations. 

Using  the  time  lag  of  the  governor,!,  the  block  diagram  for  the 
Governor  is  shown  in  Fig.  A3, 

Where:   n*       =  Rotational  speed  order 

n        =  Rotational  speed 

ner      =  n*  -n  =  Speed  error  signal 

Kf       =  Fuel  flow  rate/RPM 

T        =  Time  delay 

Wf       =  Fuel  flow  rate 

GOVERNOR  =  Kf/(l+xS) 

Then  the  propulsion  plant  is  shown  in  Fig.  A4. 
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Figure  A2.   Block  Diagram  of  Ship  Propulsion  System 
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Figure  A3.   Block  Diagram  for  the  Governor. 
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Figure  A4.   The  Propulsion  Plant  with  Governor 
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III.  HULL  EQUATION  AND  SIMULATION 


WAKE  FRACTION:   W 


The  difference  between  the  ship  speed  and  the  speed  of  advance  of 
the  propeller  Vp  is  called  Wake  speed. 
The  Wake  fraction  is  defined  as 

W  =  (V-Vp)/V 
thus      Vp=V(l-W) 
THRUST  COEFFICIENT:  Ct 
TORQUE  COEFFICIENT:  Cq 

For  the  propeller  in  open  water 


ct  =  — — rrr  (1) 

pLI       \  *  p      Til      jj      J 


Cq   =         3         \  (2) 

pD    (Vp  +n  D   ) 


Where  T  =  Thrust 
Q  =  Torque 

p  =  Water  mass  density 
Ct  =  Thrust  coefficient 
Cq  =  Torque  coefficient 
D  =  Propeller  diameter 


From  (1)  and  (2) 

2  2   2  2 
T  =  Ct.pD  (Vp+nD) 

3  2   2  2 
Q  -  Cq.pD  (Vp  +n  D  ) 


19 


The  experimental  data  on  the  Wake  fraction  versus  the  speed  of  the 
ship  is  shown  in  Fig.  Bl. 

The  Thrust  coefficient  Ct  and  Torque  coefficient  Cq  versus  the 
Second  modified  advance  coefficient  are  shown  in  Fig.  B2  and  Fig.  B3. 

Where  Vp  =  V(l-W) 

And   a  =  (n.D)/Vp2+n2D2) 

Total  ship  Resistance,  Rt,  is  a  computer  look-up  table  giving  Rt 
versus  Ship  speed  V  (fig.  B4) . 
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Figure   B2.   Thrusting  Coefficient  Ct  versus  Second 
Modified  Advance  Coefficient  o 
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IV.   COMBINED  PROPULSION  PLANT  AND  HULL  EQUATION 

A.   POWER  TURBINE  WITHOUT  GOVERNOR: 

Nonlinear  functions  such  as  Wake  fraction,  Thrust  deduction  fraction, 
Thrust  coefficient,  torque  coefficient,  ship  resistance  and  engine  torque 
were  stored  in  look-up  tables  during  the  computer  solution. 

Using  the  DSL/360  language  these  equations  were  solved  by  means  of 
an  integration  and  updating  from  table  look-up  every  0.1  sec  until  the 
steady  state  was  reached. 

Fig.  B5.  shows  propeller  shaft  speed  (rpm)  versus  time 

Fig.  B6.  shows  ship  speed  V  versus  time 

Fig.  B7.  shows  propeller  action  torque  versus  time 

Fig.  B8.  shows  shaft  torque  Qe  versus  time 

Fig.  B9.  shows  Cq  versus  time 

Fig.  BIO. shows  Ct  versus  time 
For  a  step  change  in  fuel  flow  rate 

Wf  =  7000.  +  3000.  step  (0.0) 

The  shaft  torque  decreases  and  goes  to  steady  state.   The  propeller 
action  torque  increases,  overshoots  at  10  sec.  then  goes  to  steady  state. 
The  propeller  shaft  speed  increases,  goes  to  steady  state  in  about  80  sec. 
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Figure  BIO.   Torque  Coefficient  Cq  versus  Time 
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B.   POWER  TURBINE  GOVERNING: 

Turbine  speeds  with  governor  are  shown  in  Fig.  A4  and  compared  with 
the  responses  that  occur  without  a  governor. 

For  small  step  changes  in  power  demand  at  7000  Lb/hr  fuel  flow  rate, 
the  governor,  which  has  a  rate  of  40  Lb/hr.  RPM  and  time  lag  1/^=1/5=0.2 
sec,  causes  the  power  turbine  response  to  be  reduced  from  80  sec.  to 
20  sec.  with  an  overshoot  of  3.92%  of  the  final  change  in  speed. 

RPM  command 

n*  =  120. +A.    STEP(0.) 
Now  changing 

A  =  45.  RPM  (Fig.  CI) 

A  =  40.  RPM  (Fig.  C2) 

A  =  30.  RPM  (Fig.  C3) 

A  =  20.  RPM  (Fig.  C4) 

A  =  10.  RPM  (Fig.  C5) 

Percent  Overshoot  as  a  function  of  rudder  angle  is  shown  on  Table  1. 
The  block  diagram  for  the  combined  propulsion  plant  and  hull 
equations  were  shown  in  Appendix  B. 
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V.   EQUATIONS  OF  MOTION  OF  THE  HULL 


A.   NEWTON'S  LAW  OF  MOTION: 

Newton's  Law  of  motion  for  a  rigid  body  can  be  written  as  two 
equations: 

FORCE  EQUATION 


F  =  Force  on  body  =  d (Momentum) /dt 


-  d(mUG) 

Where  U_  Velocity  of  body 


MOMENT  EQUATION 

M  =  Moment  acting  on  a  body 

d (Angular  momentum) 
dt 


d(Ift) 
dt 

Where  I  =  Moment  of  inertia 
•Sc  -   Angular  velocity 


B.   SIX  DEGREES  OF  FREEDOM: 


A 


*% 


>)e 


it 


&*> 


** 


y 
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Redrawing: 


i 


0 


z 

Y 

X 

e 


Vertical  axis,  positive  downward 
Tranverse  axis,  positive  staboard 
Longitudinal  axis,  positive  forward 
Roll  angle 


Pitch  angle 
Yaw  angle 


Positive  rotations  are  indicated 
in  the  sketch  above 


F  = 


d(mU) 
dt 


=  mU  +  X(mU) 


M  = 


_  d(H)   _ 


dt 


H+ilxH 


m:  Mass  of  ship 


U  =  ui+vj+wk    (3) 


u:  Rate  of  surging 
v:  Rate  of  swaying 
w:   Rate  of  heaving 


AO 


R=pi+qj+rk       (4) 

R:  Angular  velocity 
p  =  4>  =  Rate  of  roll 
q  =  9   =  Rate  of  yaw 
r  =  T  =  Rate  of  pitch 

~F  =  Xi  +  j  +  Zk  (5) 

X:  Hydrodynamic  force  along  X  axis 
Y:  Hydrodynamic  force  along  Y  axis 
Z:   Hydrodynamic  force  along  Z  axis 

M  =  Moment  vector  acting  on  body 
=  Ki  +  Mj  +  Nk  (6) 

K:   Rolling  moment  about  X  axis 
M:   Pitching  moment  about  Y  axis 
N:   Yawing  moment  about  Z  axis. 

Substituting  equations  (3), (4), (5),  and  (6)  into  equation  (2) 
gives 

X  =  m(u  +  qw  -  rv  -  XG(q  +  r  )  +  YG(pq  -  r)  +  ZG(pr  +  q)) 

^   Y  =  m(v  +  ru  -  pw  -  YG(p  +  r  )  +  ZG(qr  -  p)  +  XQ(qp   +  r)) 

Z  =  m(w  +  pv  -  qu  -  ZQ(p2  +   q2)  +  X^rp  -  q)  +  YG(rq  +  p))' 

K  =  Ixp  +  (Iz  -  Iy)  qr  -  m(YG(w  +  pv  -qu)  -  Z^v  +   ru  -  pw) ) 

M  =  Iyq  +  (Ix  -  Iz)  rp  -  m(ZG(u  +  qw  -  rv)  -  XQ(^   +  pv  -  qu) ) 

N  =  Izr  +  (Iy  -  Ix)  pq  -  m(XG(v  +  ru  -  pw)  -  YG(u  +  qw  -  rv) ) 

For  the  surface  ship  in  calm  water,  roll,  pitching  and  heaving  are 
taken  as  zero.   That  means  p=q=w=0  and  most  ships  have  Y  =0. 
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Then  equation  (7)  becomes: 
X  =  m(u-rv-XGr2) 


Y  =  m(v+ru+X„r  )  (8) 


N  =  Izr-m(XG(vhir)) 


C.   LINEARIZATION  OF  THE  EQUATIONS  OF  MOTION: 
1.   X  Equation 

The  nonlinear  TAYLOR  expansion  for  the  X  equation  is 


*i 


(5) 


Note: 


1.  The  TAYLOR  expansion  of  the  equation  above  includes  terms  up 
to  third  order,  since  terms  higher  than  third  order  are  not  significant. 

2.  Since  X  is  an  EVEN  function  of  v,  of   ,  of 

X(v)  =  a*  !/*+  Cln  VH+  d6  V6+  .  .  . 
X(r)  =  bi  *t  &ilH+*6n?+  '  '  ' 
X(i)  =  C,  A  CHb%  C6f+-     ' 

Then  the  odd  powers  of  v,r,<T  are  zero  i.e.  all  X  v,X  r,X/<f  are  equal 
zero. 
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Thus  equation  (5)  becomes: 

Some  of  these  terms  are  equal  to  zero  (see  Table  2) 
then  X=X =0 

2.   Y  Equation 

The  nonlinear  TAYLOR  expansion  for  the  Y  equation  is: 

y=  y+fkdK+tfv+y^+tfJ] 


+-l 


'KYV 


•J 


(7; 


NOTE: 

1.  The  terms  higher  than  third  order  are  not  significant. 

2.  Y  is  an  odd  function  of  v,r,S 

Y(v)  =  di  &  -f  ^3  ^+  As  V*+  '  '  ' 
Y(r)  =  e,  *  +  Z>>  «?+  %  ^+  '  *  * 
Y(  )  =  J,i  -f  f3/  +  £**+  •  -  * 
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Then  equation  (7)  becomes: 


+4-  yv«  **+*  a*«  «^*+  w  "*/ 


Where  the  indicated  terms  are  equal  to  zero  (see  Table  3) 


(8) 

3.   N  Equation 

The  nonlinear  TAYLOR  expansion   for   the  N  equation   is: 
HzH°  +  /Mm  +  VKV  -r/VT^  +  M<f  + 

4  3  NuufM^  t3A/m(A^t  •  -■  +$Ntun>Sl 
i6  AU^-»*+ -  '      f3N 

NOTE: 

1.  The  terms  higher  than  third  order  are  not  significant. 

2.  N  is  an  odd  function  of  v,r,  <f 

N(v)  -   #■*>+  ^^  +  05-^  '  '   * 
N(r)  =   X*  +  ^^4*5^+  '  *  * 
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Then  equation  (9)  becomes: 


The  indicated  terms  are  equal  to  zero  (see  Table  4) 
Then 

(10) 


Combining  the  third  order  TAYLOR  expansion  for  X,Y,  and  N  with  the 
dynamic  response  terms  of  the  X,y,  and  N  equation. 

X  equation    X  =  0 
Y  equation 

r 

(m-Yv)v+(mXG-Yp  r  =  Y(u,v,r,J) 

N  equation  (1) 

(mXG-N^)v+(Iz-Npr  =  N(u,v,r,<f) 
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Using  CRAMER'S  rule  to  solve  these  equations 


V  = 


Y(u,v,r,«f) 


N(u,v,r,£) 


mX  -Y 
G  r 


I  -» 

Z  r 


m-Y< 


mX  -N- 
G  v 


mX  -Y- 
G  r 


Z  r 


m-Y  • 
v 


Y(u,v,r,  ) 


mX  -N-       N(u,v,r,  ) 

\3  V 


m-Y •         mX  -Y • 
v  G  r 


mX  -N-       I  -N- 
G  v        Z  r 
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TABLE  2.   ASSESSMENT  OF  THE  COEFFICIENTS  IN 
THE  X  EQUATION 


Variable 

Coefficient 

Series   60/Model   5.1.1 

u 

AM 

Au2 
2 

V 

2 
r 

>«- 

1  V  A  U 

2 
r  au 

,2 

&  &VL 

.-  rr 
v<f    , 

vr^u 

V<5au 

m-Xu 

Xu 

1/2  Xuu 

1/6  Xuuu 

1/2  Xw 

\/<L  Xrr  +  mXj. 

1/2  Xtf 

1/2  Xwu 

1/2  Xrru 

1/2  Xtf  u 

Xvr+m 

Xv* 

Xr«T 

Xvru 

Xv  u 

Xr  u 

Xo 

0.1750 

- 

* 

Note:   No  entry  in  these  columns  means  the  coefficient  was  ignored 


From  PRINCIPLES  OF  NAVAL  ARCHITECTURE,  SNAME,  1968,  page  548, 
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TABLE  3.   ASSEMSSEMNT  OF  THE  COEFFICIENTS  IN 
THE  Y  EQUATION 


Variable 

coefficient 

Series   60/5 

,1. 

• 

V 

m-Yv 

0.509 

Y9 

mX   -Yf 

Y8 
Yl 
Y2 

V 

v5 

Yv 

1/6  Ywv 

-0.260 
-2.150 

2 

vr 

1/2  Yvrr 

-1.180 

Y3 

v*2 

1/2  YvSS 

VAU 

.   2 

VAU 

Yvu 

l/2Yvuu 

(Yr-m) 

r 

-0.0781 

YY4 

r5 

1/6  Yrrr 

-0.0461 

Y5 

2 

rv 

1/2  Yrvv 

-0.0994 

Y6 

1/2  YtSS 

rAu 

Yru 

ran* 

6 

<Jv2 

Yruu 
1/6  XT 

1/2  Y/rr 
Y«fu 

Y7 

0,Q5Q 

Sau 

Sam2 

1/2  Y*uu 

Yvr<T 

Yo 

vvS 

Yo 

0.00016 

AtC 

yyS 

(Avrf 

Y° 

1 

*No  entry  in  these  columns  means  the  coefficient  was  ignored. 
*From  PRINCIPLES  OF  NAVAL  ARCHITECTURE,  SNAME,  1968,  page  548. 
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TABLE  4.   ASSESSMENT  OF  THE  COEFFICIENT  IN 
THE  N  EQUATION. 


Variable 


coefficient 

mXG- 

■Nv 

Iz-Nr 
Nv 

1/6 

Nvw 

1/2 

Nvrr 

1/2  Nv<# 

Nvu 

1/2 

Nvuu 

(Nr- 

-mXG) 

1/6 

Nrrr 

1/2 

Nrw 

1/2 

Nr# 

Nru 

1/2 

Nruu 

N/ 

1/6  Ni^ 

1/2  N^vv 

1/2 

N<*rr 

N£u 

1/2 

N<Tuu 

Nvr£ 

No 

NS 

Series   60/Model   5.1.1 


v 
r 
v 

V5 

vr2 
v/2 

VdU 

r 

VAU* 

r 

r5 
rv2 
r*2 

S,2 
Sr2 
Jam 

vr$ 


0.0202 
•0.0750 
•0.3850 
•0.3060 


•0.0569 
■0.1010 
•1.4200 


-0.0240 


N 


MH 


-0.0003 


N9 
N8 
Nl 
N2 
N3 


NN4 
N5 
N6 


N7 


No 


*No   entry   in   these   columns  means    the   coefficient  was   ignored. 
*From  PRINCIPLES   OF  NAVAL  ARCHITECTURE,    SNAME,    1968,    page  549. 
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D.   RELATIONS  BETWEEN  EIXED  AXES  AND  MOVING  AXES: 

So  far  the  reference  axes  X,Y,  and  Z  are  fixed  to  the  moving  ship. 
In  order  to  determine  the  path  of  the  ship,  it  is  convenient  to  take  the 
fixed  axes  to  the  earth  X  ,Y  . 

In  figure  Dl.,  the  ship  is  at  the  position  A. 

The  YAW  ANGLE  y  is  the  angle  between  the  U  axis  and  the  tangent  line 
with  the  path  at  point  A. 

Then  the  following  relationships  between  the  fixed  and  the  moving 
axes  can  be  written: 


XQ(t)  =  U(t)  Cos^(t)-V(t)Sinf  (t) 


YQ(t)  =  U(t)  Sin  V(t)+V(t)Cos«/'(t) 


X0 


Figure  Dl.   Relations  Between  Fixed  Axes  in  the 
Earth  and  the  Ship  Axes. 
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E.   RUDDER  DEFLECTION: 

Let  S -   Rudder  angle 

DRATE  =  Rudder  deflection  rate. 

TDMAX  =  Time  at  maximum  rudder  deflection 


Then  the  rudder  deflection  in  turning  maneuver  can  be  shown  in 
Figure  D2 :      O, 


£   vwax 


1  ag  "kvnax 

Figure  D2.   Rudder  Deflection  versus  Time, 


Let  Rudder  deflection  =  6 


DRATE  =  Rudder  deflection  rate  =  D   /TDMAX 

max 

Then  the  function  of  the  rudder  deflection  can  be  written: 


S  =  DRATE  iRAMPCt.   )-RAMP(t-t    )1  -r=r 
lag  max  J  57. 
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F.   THE  TURNING  PATH  OF  A  SHIP: 

There  are  4  phases  of  a  turn: 


TACTICAL     P/AME7»? 

lit  Chase  j  \ 

!I         [<']     EXECUTE    POINT 
API?-  PHASE 

I 

Figure  D3.   The  Turning  Path  of  a  Ship. 

Approach  phase  V=0,  r=0,v=0,  r=0 
First  phase  of  turn  vfO,  r  ^  0,  v  =  0,  r=0 
Second  phase  of  turn  v=^0,  r  ^  0,  v  ^  0,  r/0 
Third  phase  of  turn  v=0,  r=0,  v^O,  r^O 
Using  dimentional  equation  (1)  and  introducing  the  rudder  deflection 

30  .   The  characteristics  of  transient  phases  of  a  turn  are  shown  in 

Figures  D6,  D7,  D8,  and  D9 . 


A-B  Second  phase 
B-C  Third  phase 


V7*0,  v  ±  0,   r^O,  r^O 
v^O,  v  =  0,  r^O,  r^O 


Now  using  different  rudder  deflections  =35  ,30   the  turning  path 
of  a  ship  is  shown  in  Figures  D4  and  D5.   When  using  the  bigger  deflection 
the  steady  turning  radius  is  smaller. 
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Figure  D6 •   Angular  Acceleration  about  Z  Axis. 
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Figure     D7-      Angular  Velocity  about   Z  Axis 
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Figure   D8.   Angular  Acceleration  about  Y  Axis 
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Figure  D9.    Angular  Velocity  about  Y  Axis. 


25.  r. 


58 


VI.   SPEED  CONTROL  IN  TURNS 

A.   DYNAMIC  BEHAVIOR  OF  THE  PROPULSION  PLANT: 

The  propulsion  plant  and  turning  path  were  observed.   Now  combine 
both  of  them  and  add  feed-back  of  ship  speed  with  Gain  Kl  in  order  to 
keep  the  ship  speed  constant  in  turns. 

The  block  diagram  for  this  part  can  be  drawn: 


V*   _  Vex 

— sO 


K± 


n* 

'■?  Wei 

s 

WO  PULSION 
PL  AWT 

T> 

HULL 

VI 

v-6-  v 


Figure  El.   Feedback  Control  Velocity. 

Where:   V*:  Reference  ship  speed. 

V  :  Ship  speed. 

n*:  Shaft  angular  speed  command. 

U  :  Velocity  about  X  axis. 

VA:  Velocity  about  Y  axis 

Kl:  Ship  speed  gain 

n  :  Shaft  angular  speed 

6  :  Rudder  angle 

The  block  diagram  of  this  system  was  shown  in  Appendix  B, 
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B.  SIMULATION  METHODS 

1.  Determine  all  initial  conditions  in  the  steady  state 

RUN  1. 

Let  propeller  shaft  speed  command  n*  be  constant 

n*  =  130  RPM 
And  no  feedback  of  ship  speed 

Kl  =  0 
Then  determine  the  steady  state  speeds  (i.e.,  determine 
the  component  of  initial  condition  ship  speed  on  X  axis 
and  on  Y  axis  (UIC,  VIC). 

RUN  2. 

Now  the  propeller  shaft  speed  command  n*rl30  rpm  making 
Kl  =  8,16,24,32  respectively.   Then  determine  the  steady 
state  in  each  case  (i.e.,  find  all  IC's  in  the  steady 
state) .   This  guarantees  that  the  loop  is  in  steady 
state,  since  the  initial  conditions  on  other  variables 
might  still  cause  a  transient. 

2.  Put  in  a  step  change  in  the  ship  speed  with  2  kts.  increase  and 
delay  time  10  sec.   Then  V*  becomes: 

V*  =  20.  +  2.0  step  (10.) 

C.  SPEED  CONTROL  STUDIES 

Using  the  inital  conditions  of  RUN  1  and  2,  a  rudder  command  was 
applied  at  t  =  0,  causing  a  speed  decrease  due  to  added  drag.   At 
t  =  10  sec,  the  commanded  speed,  V*  was  increased  by  2  knots,  stepwise. 
This  caused  a  speed  increase.   The  effects  of  changing  the  loop  gain, 
Kl,  is  shown  in  Figure  E2. 
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D.   DISCUSS  THE  CURVES 

1.  Ship  resistance  Rt,  Wake  fraction  W,  Thrust  fraction  t 

The  ship  speed  change  was  small,  so  that: 

a.  The  ship  resistance  Rt  =  K.V   (see  Fig.  B4)   K  =  constant 
then  the  form  of  Rt  was  the  same  as  V  (Fig.  E3  curve  1) 

b.  The  Wake  fraction  was  constant  (Fig.  Bl) 
then  W  =  (Fig.  E3  curve  2) 

c.  In  this  interval  the  slope  of  the  thrust  fraction  was 
negative.   Then  when  V  decreases  to  increases  and  as  V 
increases  t  decreases  (Fig.  E3  curve  3) 

2.  a  curve,  torque  coefficient  Cq,  thrusting  coefficient 

VP  =  (1-W)  V 
=  (1-.045)  V 

=  -655V 

,  nD 

and     a     = 


„  2   22 
Vp  -n  D 


The  a  curve  is  shown  in  Fig.  E4  curve  1 

a  changed  from.^to  .520  then  the  slope  of  Ct  and  Cq  curve  was  constant, 
Then  the  Ct  and  Cq  curve  have  the  same  form  as  the  <T  curve. 

3.   Propeller  action  torque  Qp,  shaft  torque  Qe  and  propeller  angular 

speed  n 

3   2   2  2 
Qp  =  Cq.p.D  (VP  +n  D  ) 

The  Qp  curve  was  shown  in  Fig.  E5  curve  2. 

The  Wf  curve  and  n  curve  are  shown  in  Figures  E6  and  E7. 
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Figure  E4.   Second  Modified  Advance  Coefficient    ,  Torque  Coefficient 
Cq,  Thrust  Coefficient  Ct  versus  Time. 
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Figure  E5.   Propeller  Action  Torque  Qp  and  Shaft  Torque  versus  Time. 
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Figure  E6.   Fuel  Flow  Rate  Wf  versus  Time. 
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VII.  CONCLUSION 

This  thesis  has  presented  the  combined  model  of  a  propulsion  plant 
and  hull.   DSL/360  LANGUAGE  was  found  to  be  powerful  for  simulating  ship 
motion. 

In  turbine  powered  ships,  it  is  necessary  to  predict  the  propeller 
response  to  changes  in  power  demand. 

This  problem  has  been  investigated  but  the  studies  have  been  con- 
centrated on  the  effect  of  dynamics  of  the  ship.   This  work  is  concerned 
both  with  the  dynamic  of  the  propulsion  plant  and  ship  dynamics  of  the 
ship. 

It  has  been  shown  that  the  complex,  nonlinear  dynamics  of  propulsion 
plant  and  hull  can  be  adequately  represented  by  an  all  digital  computer 
model.   The  model  permits  studies  of  the  internal  dynamics  of  the  pro- 
pulsion plant,  as  well  as  the  external  dynamics  of  the  hull. 

It  has  been  shown  that  the  model  can  be  incorporated  in  a  feedback 
control  system  and  its  effects  included  in  studies  of  control  system 
performance. 
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VIII.   RECOMMENDATIONS 

The  combination  of  propulsion  plant  and  hull  in  a  single  computer 
model  provides  a  realistic  tool  for  many  studies.   The  following  are 
topics  suitable  for  future  investigation. 

a.  Course  keeping,  station  keeping  and  replenishment  at  sea. 

b.  Effect  of  propulsion  plant  dynamics  on  maneuvering  in  sea 
states  (or  regular  waves) . 

c.  Effect  of  sea  state  conditions  on  the  operation  of  the 
propulsion  plant. 

d.  Discuss  variation  in  steady  state  speed  with  increase  command 
as  shown  in  Fig.  CI  -  C5.   Also  discuss  effect  of  governor 
loop  on  steady  state  speed,  Fig.  E2. 
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APPENDIX  A 


Table  Look-up  for  Ship  Resistance  versus  Speed; 


300     - 


250   - 


FUNCTION  TABR(V) 

DIMENSION  VT(13),RTT(13) 

DATA  VT/0. 0,2. 5, 5. 0,7. 5, 10., 12. 5, 15., 17. 5, 20., 22. 5, 25. 5, 127. 5, 30./ 

DATA  RTT/2*0.0, 7000. ,13000. ,25000. , 39000. ,57000. ,80000. ,1103000. , 

136000., 173000., 225000. ,280000./ 

IF(V.GT.30)  GO  TO  1 

DELV  2.5 

N  IFIX(V/DELV)  1 

SLOPR (RTT (N  1) -RTT (N) ) / ( VT (N-l) -VT (N) ) 

TABR  SLOPR* (V-VT(N)) -RTT (N) 

RETURN 

1  TABR  280000. 

RETURN 

END 
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